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Phonon-bottleneck enhanced magnetic hysteresis in a molecular paddle wheel 

complex of Ru2 + 
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The ruthenium based molecular magnet [Ru2(D(3,5-Cl2Ph)F)4Cl(0.5H2O)-C6Hi4] (hereafter RU2) 
behaves as a two-level system at sufficiently low temperatures. The authors performed spin detec- 
tion by means of single-crystal measurements and obtained magnetic hysteresis loops around zero 
bias as a function of field sweeping rate. Compared to other molecular systems, Ru2 presents an 
enhanced irreversibility as shown by "valleys" of negative differential susceptibility in the hystere- 
sis curves. Simulations based on phonon bottleneck model are in good qualitative agreement and 
suggest an abrupt spin reversal combined with insufficient thermal coupling between sample and 
cryostat phonon bath. 



Molecular magnets have been explored intensively in 
the recent years for their potential application in informa- 
tion technology and as test beds for quantum mechanics 
at macroscopic scale. These systems are tunable, iden- 
tical, two- or multilevel systems that can be produced 
in large numbers and therefore are potential candidates 
for qubit implementation in quantum computing algo- 
rithms 1]. In crystalline form, the molecules are rela- 
tively well isolated from each other and thus the quan- 
tum physics of their spin is not strongly affected by the 
collective nature of the sample. Large spin molecular 
magnets have demonstrated key quantum effects, such 
as quantum tunneling of the magnetization through an 
anisotropy barrier 0, [H and Berry phase interference [H • 
Low spin molecules have shown interesting effects when 
tuning phonon dissipation [j| that provide increased re- 
laxation and decoherence times for the molecular qubits. 

In the case of low-spin molecular systems, the 
anisotropy barrier against spin reversal is often very 
small. However, magnetic hysteresis has been recently 
measured for a series of low-spin molecules (Vj.5 S=l/2, 
Feio S=l, V 6 S=l/2) 0, @, 0, i and explained in the 
frame of the phonon-bottleneck (PB) model [HQ. In 
this model, the phonons' low heat capacity keeps the lat- 
tice and the spins at the same temperature T s , slightly 
different from cryostat temperature T, as a result of in- 
herent limited sample thermalization in the experimental 
setup. Cycling of an external magnetic field induces a de- 
lay (hysteresis) in the dynamics of the spin temperature 
which is observed by measuring the magnetic moment of 
a single crystal. Other theoretical interpretations treat 
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the lattice at equilibrium with the cryostat at a tem- 
perature T ^ T s |10j. Though magnetic hysteresis can 
be simulated for a given bath temperature, these models 
cannot explain the observed dependence on sample ther- 
malization [5[ and the picture of phonons having very low 
heat capacity at low temperatures Q . On the other hand, 
the PB model presented in Ref. @ gives a good quantita- 
tive agreement with the experiments. The present study 
introduces a particular aspect of the phonon-bottleneck 
phenomena as observed in the Ru2 paddle wheel molec- 
ular system (structure in Fig. [I]). Single crystal mea- 
surements reveal enhanced magnetic irreversibility with 
"valleys" of negative differential susceptibility followed 
by rapid increase of magnetization. The observed dy- 
namics can be explained quantitatively in the frame of 
the PB model for an abrupt spin reversal in zero field, 
followed by a rapid gain in phonon energy. 

The [Ru 2 (D(3,5-Cl2Ph)F)4Cl(0.5H 2 O)-C 6 Hi4] molec- 
ular crystal was synthesized following the methods de- 
scribed in Ref. [11]. It contains a magnetic Ru^ + core 
with Ru— Ru bond length of 2.36 A, surrounded by bridg- 
ing ligands in a paddle wheel pattern. Each molecule is 
formed of a dinuclear core which allows the ground state 
of the molecule to be described as a single spin. The crys- 
tals are tetragonal, and the Ru-Ru axes of all molecules 
are parallel to one another and aligned with the c axis 
(001) of the unit cell. The paramagnetic units are sepa- 
rated by 11.2 A (nearest neighbor distance) so that the 
dipolar interaction between molecules is very weak. The 
compound can be synthesized in relatively large crystal 
sizes with low-index faces suitable for single crystal ex- 
periments. Magnetic resonance studies [Tl| indicate a 
total spin S = 3/2 with the lowest two levels S z = ±1/2 
located ~ 200 K below the S z = ±3/2 levels. Conse- 
quently, the Zeeman splitting in a longitudinal field cor- 
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FIG. 1: The [Ru 2 (D(3,5-Cl2Ph)F)4Cl(0.5H2O)-C 6 Hi4] 
molecule with nondisordered atoms shown as displacement 
ellipsoids at the 30% probability level (H atoms are omitted 
for clarity). The two Ru ions form the magnetic core of the 
molecules, synthesized in a single crystal with tetragonal sym- 
metry. The Ru-Ru axes of all molecules are parallel to one 
another and aligned with the c axis (001) of the unit cell. 
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FIG. 2: Experimental (a) and simulated (b) hysteresis curves 
for a magnetic field along the 001 axis of the sample measured 
at various temperatures and magnetic field sweeping rates. 



FIG. 3: Hysteresis curves for a magnetic field along the 110 
axis of the sample: (a) experimental data measured at various 
temperatures and magnetic field sweeping rates; (b) simulated 
cycles in the same conditions as in (a). 



responds to a spin 1/2 system, whereas in a transverse 
field it is twice as large due to increased mixing with the 
upper levels (Tlj. 

The effective Hamiltonian of a two-level system can 
be written as H e ff = —ea z /2 — Ao<Xe/2, where a ZlX are 
the Pauli spin matrices, An is the level repulsion, and 
e = g/ifii? is the spin bias (/is is the Bohr magneton 
and H is the external magnetic field, always || to z-axis 
of H e ff). The above mentioned effect of a transverse 
field on the Zeeman energy can be described by a gyro- 
magnetic factor g = 2 when H || 001 and g = 4 when 
H || 110. Thus, the energy separation can be written as 
Ah = \/ Aq + e 2 , which is also the energy Uu> = Ah of 
phonons responsible for spin thermalization. 

To analyze the role of phonon bottleneck on the spin 
dynamics of the Ru2 molecule, we performed single- 
crystal magnetic measurements using a micron-sized su- 
perconducting quantum interference device [3| at very 
low temperatures (^1 K) when the molecule is purely a 
two-level system. Hysteresis cycles in field swept between 
±1.4 T at constant rate are presented in Figs.[2k[3K (f° r 
clarity, only one half of the cycle is shown; the other half 
being symmetric against the origin). As in the case of 
other low-spin molecules [(| 0, § , an irreversibility is 
occurring even though no anisotropy barrier against spin 
reversal was detected. 
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In the case of the R112 molecule we notice a large 
spin flip during passage through the origin, as shown by 
an abrupt ascending branch that deviates from the re- 
versible behavior and can reach ~90% of the total mag- 
netization (see data for T = 40 mK, r = 8 and 2 mT/s in 
Figs.EK and [3^ respectively). This is the signature of a 
small level repulsion Ao, which in the absence of thermal 
broadening and for sufficiently small spin temperatures 
is observed as an abrupt spin flip (see also Ref. 0). We 
note that the nearly full spin reversal is still done in a 
quantum adiabatic fashion because of fairly low sweeping 
rate (no Landau-Zener excitations). 

When in thermal equilibrium and for Ao large enough 
to prevent nonadiabatic excitations, the molecule spins 
follow H adiabatically through zero field showing a spin 
up - spin down transition, as observed at low sweep rates 
(see Fig. [2] for T=0.4 K and r = 4mT/s). However, 
for sufficiently large sweeping rates the number of avail- 
able phonons at energy Ah{H) can decrease if sample 
thermalization is not perfect. In Fig. [3] a small effect 
is visible already at T=0.4 K and r = 8mT/s. At low 
temperatures phonons have much less heat capacity than 
spins 9] and the spin-phonon coupled system can be 
treated as being at the same temperature T s , defined 
by n\/n2 = exp(A# /)cbT s ), with 77,1,2 the out of the 
equilibrium occupation numbers of the ground and ex- 
cited spin levels, respectively. Our experiments show a 
steep decreasing field branch [see arrows for the first four 
curves in Figs.^k and[3k] indicating T s < T, where T is 
the cryostat temperature. This behavior is due to a lack 
of thermal coupling between the lattice and the phonon 
reservoir of the cryostat leading to an insufficient num- 
ber of available phonons for spin thermalization. This 
behavior persists after passing through zero field as well. 
Quantitatively, this phenomenon can be described by the 
PB model @, Q in which deviations from equilibrium are 
given by xt h = l/x—1 with x = (ni - n 2 )/{n leq - n 2eq ) 
and spin relaxation time @ : 

tanh 2 (A ff /2fc B T) 
T H = a -~ 2 ■ (!) 

H 

The parameter NT p h/3Auj is proportional to 

molecule density N, phonon velocity v and phonon-bath 
relaxation time t v % (depending on the sample thermal- 
ization to the external bath) and inversely proportional 
to the spectral width of the levels Au>. 

In contrast to other molecular systems, the Ru2 spin- 
flip induces a large increase of phonon energy that is vis- 
ible through a valley of negative differential susceptibil- 
ity. The magnitude of the effect and its dynamics depend 
on the competition between two field-dependent param- 
eters, the relaxation rate th and the time Ah /v spent at 



a given two-level energy separation (with v = dAn / dH) . 
The PB model can describe accurately this phenomenon 
as one can see from the simulations of Figs. Wp and [3Jd. 
Following the abrupt spin reversal, the spin system is 
"scanning" again the phonon energy distribution and 
available resonant phonons are used to excite the spins. 
Subsequent spin deexcitations will generate phonons with 
slightly higher energy since Ah has been increased mean- 
while. In a short time, these dynamics increase the 
phonon energy density which compromises the magneti- 
zation rapidly and stimulates emission of further phonons 
@ . The gain in phonon energy, combined with the imper- 
fect thermal coupling of the sample, leads to an increased 
spin temperature. Thus, in this field range the phonon 
bottleneck acts in the opposite way, that is, the phonons 
cannot "escape" from the sample into the external bath. 

The simulations using the PB model require two fit 
parameters, namely, the zero-field level repulsion Ao 
and the sample thermalization parameter a. We find 
A = 15 mK and a = 1 and 3 sK 2 for H || 001 and 
H || 110, respectively. The measurements of Figs. [2^ and 
13^, were done in two separate runs with inherently differ- 
ent sample thermalizations to the cryostat. The resulting 
different values of a indicate an intermediate thermaliza- 
tion. As a comparison, in related experiments on V15 
molecule @, a was 0.09 and 130 sK 2 for a sample well 
thermalized and well isolated, respectively, from the cryo- 
stat. 

The obtained value of Ao is several times smaller than 
the value found for V15 (A = 80 mK d). 

This causes a 

more complicated situation near zero field, where the spin 
temperature is strongly reduced. In this case, the PB 
model assumption T s < T is no longer valid, which affects 
the agreement between simulations and experiments. 

In conclusion, we introduced a particular aspect of 
the phonon-bottlencck phenomena as detected by single- 
crystal magnetic measurements of the RU2 paddle wheel 
molecular system. Ru2 is a member of a class of molecu- 
lar magnets with a very simple (diatomic) core that ex- 
hibits this phenomenon; all the previous cases involved 
much more complex spin metal clusters. Enhanced mag- 
netic irreversibility with "valleys" of negative differential 
susceptibility followed by rapid increase of magnetization 
is observed. These dynamics are due to an abrupt spin 
reversal in zero field, followed by an increase in the en- 
ergy of phonons, and can be explained quantitatively in 
the frame of the PB model. 
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